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An analytical solution for light field modes in waveguides with nonideal cladding 
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We have obtained an analytical solution for the dispersion relation of the light held modes in the 
nanowaveguide structure. The solution has been analyzed for the planar waveguide with metama¬ 
terial claddings and dielectric core. The analytical solution is valid within the broadband spectral 
range and is confirmed by existing numerical calculations. The developed theoretical approach opens 
vast possibilities for the analytical investigations of the light fields in the various waveguides. 
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Introduction. Control, manipulation and spatial con¬ 
finement of light fields are the main objectives of photon¬ 
ics where various waveguides and resonator structures 
play a major role [D-Q. The growing interest in this 
field of research is spurred by the invention of diverse 
nanowaveguide systems in optics, photonics, plasmonics 
and its applications to manipulate the single photon fields 
M- Great expectations are associated here with using 
new materials, such as photonic crystals, metamateri¬ 
als, graphene, and others [fi [§-16 causing the detection 
of entirely new properties of light. The strong stream 
of these experimental investigations stimulates elaborat¬ 
ing the new theoretical approaches which could provide 
a clearer physical picture and high mathematical preci¬ 
sion for the description of the light fields in the new and 
well-known waveguides. 

Theoretical study of the light field localized in the 
waveguide structures is a long-term theoretical prob¬ 
lem a Work on the spatial confinement of 

the free light fields into the waveguide with the finite 
cross-section shows a dramatic change in the fundamen¬ 
tal properties of light. The transversely confined pho¬ 
ton wave packet propagating along z —direction acquires 
the spectral dispersion in the form of Hamiltonian func¬ 
tion huj = \J ( hck z ) 2 + W n corresponding to the rela¬ 
tivistic particle with a finite mass too j23l425|. where 
W n = ( to 0 c 2 ) 2 = (hckn) , w = 2ttc/X , c is a speed 
of light, A its wavelength in the vacuum, h - Planck 
constant, /c~-longitudinal wave-number of the n-th light 
mode, and k^ is a transverse wave number determined 
by the specific properties of the waveguide (see later). It 
is worth to note the recent experiment [2(| where effec¬ 
tive nonzero photon mass revealed itself in a possibility 
of Bose-Einstein condensation of photons in the optical 
microcavity. 

In practice, the waveguides have various geometries, 
cores and materials of the waveguide claddings 21] which 
are characterized by different spectral properties of per¬ 
mittivity e(u>) and magnetic permeability n(ui). Geomet¬ 
rical and physical properties of the waveguides determine 
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FIG. 1. (Color online) The three-layer planar waveguide 
consists of the core layer of transverse dimension L = a + b 
, and e co (w), p co (w). The thicker claddings are located at 
x = b and x = —a and made of e±(w),/x±(w). The light 
modes propagate along the waveguide axis 2 . 


highly specific behavior of the light modes il> n at the in¬ 
terface between the claddings and cores identifying the ef¬ 
fective cross-section of the confined light modes, the con¬ 
crete mass to 0 of photon and its minimum frequency ui n , 
respectively. In particularly, co n = and Too = for 
the simplest planar waveguide with the transverse width 
L and a perfect metal cladding. 

The conhned light modes in the arbitrary waveguides 
are described by the transcendental dispersion equations 
HEH that do not have analytical solutions even in the 


simple cases. Unfortunately, the transcendental equation 
is also impracticable for intuitive analysis of the light 
properties. In turn, the understanding of the eigen mode 
dispersion relation is required for the further study of the 
principal information about the localized electromagnetic 
fields such as the phase and group velocities, transverse 
shape and propagation length of the low-losses modes, in 
particular for the planar or nanowire arrays and for the 
nanoplasmonic, V-groove structures, etc. (see ana and 
references there). 

In this letter we present the transparent analytic solu¬ 
tion for the dispersion relation of the confined light modes 
in the nanowaveguide structure. Here, we have used the 
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mathematical method of Kacenelenbaum B. [28j and de¬ 
veloped the perturbative analysis of term W n related to 
the photon mass m 0 ^ 0 manifested in the Hamiltonian 
function of photon in any type of waveguide regardless 
of the physical features of its walls. To the best of our 
knowledge this approach is accomplished here for the first 
time to the analytical solution which can be applied with 
high accuracy to the theoretical study of the real waveg¬ 
uide structures. 

The perturbative approach. We demonstrate our ap¬ 
proach for the well-known three layer planar waveguide 
consisting of the thin slab and two much thicker claddings 
as depicted in Fig. [T] This is the simplest non-trivial sys¬ 
tem where the salient features of the current approach 
can be clearly demonstrated. The eigen confined light 
modes are expressed here through the Hertz potentials 


[ A n e~ k+X , x>b, 

i’n = < B n smk^x + C n cosk^x, (1) 

{D n e k - X , x < -a, 

where e l i k ^ z ~ ult ) i s omitted (for notaional convenience), 
k z is the wavenumber parallel to the interface, k^ is the 
transverse wavenumber in the core, k 2 + k^ 2 = kQ£ co p , co , 
k± are the transverse wavenumbers in the claddings k± = 
k 2 — k^e±p± (where indexes ”co” and ” ± ” denote the 
core and the claddings, see also Fig[l|). In this letter we 
restrict our attention to the transverse magnetic (TM) 
modes (the transverse electric (TE) modes can be studied 
in a similar way) for which the electric and magnetic 
components of the modes are written as follows: 


\ uzuj; ' 


From the boundary conditions for the fields at the two 
interfaces x = —a and x = b we can get the dispersion 
relation in the form of transcendental equation similar to 


Our aim is to develop the perturbative method fll] for 
the solution of the original wave equation A %/; n +k^ 2 ^>„ = 
0 by taking into account W n = (hckn) and the fun¬ 
damental Hamiltonian form of the dispersion relation 
lo 2 = ( ck ^) 2 + ( ck z ) 2 . Here, the transverse wavenum¬ 
ber square k^ 2 and the eigen mode ip n are decomposed 
in a series expansion with respect to a small parame¬ 
ter | q | < 1 (the physical meaning of the parameter q is 
determined by the parameters of the material claddings 
that will be discussed later): 

k n 2 = k n\o) + Q k nt 1 ) + - + ^ k n (j) + - . ^ 

4>n = Ipn' 1 + qipn' 1 + + q j tpn' > + ■■■ 

By inserting Eq.0 in the wave equation and equating 
the terms of the same order to the small parameter g, we 
obtain Ai/'i 0 ' 1 + = 0 and A + 

kn (o)^’" 1) = 0- One can (i) multiply the first equation by 
ipi, the second by ipo, (ii) integrate the equations over the 
waveguide cross section S by taking into account Green’s 
theorem and normalization k^ 2 ^ f dSi/jn^ 2 = 1, and (iii) 
subsequent subtraction of two integrals leads to: 



where C is the contour of the cross section S', A is a 
unit vector normal to the interface between the core and 
claddings. 

The solution k^ 2 ^ is expressed through the bound¬ 
ary values of the functions t/ffl and ipn^ on the con¬ 
tour C. In order to find the functiond t/>o, ipi on the 
C , we apply the Rytov-Leontovich boundary condition 
successfully used for studing the light effects on the in¬ 
terface between two different materials 31). This condi¬ 
tion determines the relation between the magnetic and 
electric fields of light mode through the impedance C± of 
claddings at the boundary contour C: 


tan (kn(a + b) — mr) = 


(& + fc) 


■co GO 




( 3 ) 
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Unfortunately, the well-known theoretical approaches 
permit only the approximate solutions of the transcen¬ 
dental equation in three limiting cases: i) near cuttoff, 
ii) at short-wave limit Lw/c 1, and iii) in the strongly 
asymmetrical case (see, for example 0 )- However, the 
most practical systems are characterized by the inter¬ 
mediate sizes Llo/ c < 1 (various nanowaveguide struc¬ 
tures etc.), by different geometries and material proper¬ 
ties where only the numerical methods are possible for 
solving the complex transcendental equations 21]. 


E n 
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( 6 ) 


where C±( w ) = t L ±i UJ )/ £ ±( UJ ) and it is assumed that £± 
has the same order of smallness as parameter q. 

Substituting Eqs. ©, 0 in Eq. (0 and equating 
the zero-order terms, we find = 0, this corre¬ 

sponds to the TM mode in the waveguide with perfect 
metal cladding. Here, we have i/jn' 1 = B n sin(fc^ 0 ^a;) + 
C n cos(k^ 0 ^x), where B n and C n are determined from 
the boundary conditions on the countour C and by nor¬ 
malization condition of ipn\ k n(o) ~ (ufa) • By per- 
forming the same calculations with Eqs. 0, 0 in Eq. 













3 


© for the first order of smallness £± and by taking into 

account that = 0 , we get for : 

o o 
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(7) 


Using the fact that the impedance is independent of the 
transverse coordinates and substituting ipn\ and i/jn' 1 in 
Eq. (©, then we find the analytical expression for the 
wave number in the asymmetric waveguide up to first 
order of the pertubation expansion: 


^ = (s) 2 -K++<-) 2 Sf. 

k 2 z = k>e d p d -(^) +(C+ + C-) 2 tf £L - 

In case of the symmetric waveguide (£+ = £_ = £ c i and 
b = a), the dispersion relation of even modes is given by: 


fci = (ir) 2 -Cc^, 

kl = kle d n d -(^) 

and for the odd field modes we have: 


k\ 


= (■?- 
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(9) 


( 10 ) 


where n = 0 , 1 , 2 ,... defines a set of waveguide modes. 

The series of the analytical solutions m-m are the 
main result of this work. It is interesting to compare 
the solutions with the numerical results obtained for the 
waveguides with realistic physical parameters. 

Comparison with the numerical results. The intensive 
numerical studies of the light field propagation in the 
planar nano waveguide with the metamaterial claddings 
have been performed in the recent works [32| . l33| . This 
analysis required quite large computing resources. In 
case of the metamaterial-dielectric interface, the light 
held can be pushed out from the metamaterial into the 
dielectric core when dielectric permittivity and magnetic 
permeability of the metamaterial are both negative. As a 
result, the low-losses held modes are excited on the dielec- 
tric/metarnaterial interface 


34(, 


j. Nanowaveguides 
[27l ] with such nretanraterial claddings seem to be promis¬ 
ing for realization of highly confined low-losses modes. 

We compare the numerical results for TM even modes 
[32 ;, [33| with our analytical solution Eq. ©. The 
metamaterial claddings were described by the Drude-like 
model of permittivity and permeability: 


e(w )/£ 0 = 1 - We/(w(w + *7e)), 
p(ui)/n 0 = 1 + Fuj 2 /(u}% - w(w + *7m)), 


( 11 ) 


where the electric 7 e = 2.73 • 10 13 s 1 and magnetic 
7 m = 7 e damping rates are much less than the carrier 



FIG. 2. (Color online) Spectral behavior of the refractive 
index \n c i\ (blue line), and the impedance |£ c i| (red dashed 
line) of the metamaterial as functions of w/w e . 


1.2 

1.0 

0.8 

0.6 

0.4 

0.2 


FIG. 3. (Color online) The behavior of Re{k z } for the TM 
modes: TM n =i, TM n= 2 depending on the frequency; the an¬ 
alytical solutions (solid curves) and the numerical (dashed 
curves) solutions are presented here for the transverse size 
2a = Anc/uJe « 275 nm, and e co = 1.3, p co = 1. 



frequency of interest; 7 e ^ m to, u> e = 1.37 • 10 16 s -1 is a 
plasma frequency of the material, w 0 = 0 . 2 w e is a binding 
frequency, and F = 0.5 is a geometrical factor accounts 
for the magnetic oscillation strength. 

The spectral properties of the impedance Cci(w) and 
the refractive index n c i(oj) ( n 2 cl {<jj) = e(u))p(oj)) for used 
metamaterial cladding are presented in Fig. [2j As it 
is seen in Fig. [2j the condition | ( c j \ < 1 is satisfied for 
the frequencies w < 0.85w e (except small area around 
to = 0.2w e ), so the perturbative expansion Eq. a is 
valid for most of the spectral range. It is seen in Figs. 
HED that each TM mode has its spectral domain where 
the good match between the analytical and the numer¬ 
ical solutions occurs. In particular, for the first TM 
mode (n=l), we have found the analytical solution fc“ n 
for the longitudinal wave number coincides with the nu¬ 
merical one k™ um with precision ( kz | —( < 10 -2 

for the spectral range w < 0.25w e (except for a small 
spectral range ~ 0.04w e around w = 0.2 ui e ) where the 
Rytov-Leontovich condition holds. In the spectral range 
0.25w e < uj < 0.85w e , the precision of the wave number 
remains high * | k nl m z \ —- < 5 • 10 ~ 2 (i.e. in 5 times 
smaller, except the small area around cut-off frequency 
0.4w e )- So as it is seen in Fig.4, the difference between the 
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FIG. 4. (Color online) The behavior of Im{k z } for the 
TM modes: TM n= i, TM n= 2 , depending on the frequency. 
The analytical solutions (solid curves) and the numerical 
(dashed curves) solutions are given for the waveguide width 
2 a = A-Kc/uJe ~ 275 nm. The enlarged scale inset shows 
the sharp dips in the attenuation coefficients close to the 
frequency ~ 0.2oj/u} e for the both light modes (the curves 
demonstrate perfect coincidence between the analytical and 
the numerical solutions for these dips). 


analytical and the numerical solutions for attenuation co¬ 
efficient is negligible within the broad spectral range of 
the figure. The similar high accuracy of the analytical 
description occurs for the attenuation coefficient of the 
second TM-mode (n=2). 

As seen in Fig. [2 the small dips in the attenuation 
coefficient for the both TM light modes demonstrate an 
emergence of the low-losses field modes. Our prelimi¬ 
nary analysis of the analytical solution Eq. m in this 
spectral range shows the spectral dip parameters (width 
and depth) are highly dependent on the waveguide trans¬ 
verse size, the field interference effects and the intensity 
distribution in the waveguide cross section. A clearer 
understanding of the main features in the suppression of 
the low loss mode attenuation requires further detailed 
studies for various parameters of the waveguide. 

We have also performed the accuracy analysis of the 
analytical solutions for the TM mode (n = 1). By tak¬ 
ing into account the spectral behavior of the impedance 
and the refractive index we conclude the following state¬ 
ment. The solution Eq. © has a large accuracy within 
the quite broad spectral range where the impedance is 
sufficiently low | Cci | < 1 and moreover it occurs even for 
the | (d | ~ 1 provided the refractive index of cladding 
materials is sufficiently large \n c i\ > \ ( c i \. Herein, the 
higher orders of the perturbation terms in series Eq. © 
make a negligible contribution in the analytical solution. 

Conclusion. In this letter we have analytically solved 
the critical problem of the dispersion relation of the field 
modes in the waveguide where light modes are strongly 
localized. Herein, we have developed the perturbative ap¬ 
proach for light propagation in the nanowaveguide struc¬ 
tures that is based on the properties of Hamiltonian func¬ 
tion of the light fields in the transversally confined waveg¬ 
uides. The perturbation approach is constructed for the 


decomposition of the transverse wavenumber in a small 
parameter q determined by the physical properties of 
the waveguide claddings. The obtained analytical results 
have demonstrated a very good accuracy compared with 
the well-known numerical results within the broadband 
spectrum. 

Thus, we have developed the theory demonstrating a 
unique predictive opportunity for the analytical study of 
the complicated waveguides characterized by small and 
even moderate impedance of their claddings. The ob¬ 
tained analytical solution is more convenient and practi¬ 
cal for studying the light field modes in the waveguides, 
with various physical and geometrical parameters, which 
are important components in the conventional optics, 
nanooptics, integrated QED- circuits as well as in the 
nanoplasmonics. 
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